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Abstract: The study of time lags between spikes in Gamma-Ray Bursts light curves in different energy 
bands as a function of redshift may lead to the detection of effects due to Quantum Gravity. We present 
an analysis of 15 Gamma-Ray Bursts with measured redshift, detected by the HETE-2 mission between 
2001 and 2006 in order to measure time lags related to astrophysical effects and search for Quantum 
Gravity signature in the framework of an extra-dimension string model. The use of photon-tagged data 
allows us to consider various energy ranges. Systematic effects due to selection and cuts are evaluated. 
No significant Quantum Gravity effect is detected from the study of the maxima of the light curves and a 
lower limit at 95% Confidence Level on the Quantum Gravity scale parameter of 3.2 x 10 GeV is set. 



Introduction 

Two approaches are followed trying to elaborate 
a Quantum Gravity theory, which would make 
the link between particle physics and cosmology: 
String Theory and Loop Quantum Gravity (which 
will not be discussed here, see [1]). In String 
Theory, gravitation is considered as a gauge in- 
teraction and Quantum Gravity effects result from 
graviton-like exchange in a background classical 
space-time. 

In this paper, we follow [2, 3] and we consider a 
model [4] based on String Theory where photons 
propagate in the vacuum which may exhibit a non- 
trivial refractive index due to its foamy structure 
on a characteristic scale approaching the Planck 
length I ~ nip\ anck . This implies a light group 
velocity variation as a function of the energy of the 
photon E: v(E) = c/n(E), where n(E) is the re- 
fractive index of the foam. Generally, the Quantum 
Gravity energy scale Eqg is considered to be close 
to the Planck scale. This allows us to represent the 
standard photon dispersion relation with E/Eqq 
expansion: 
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where £ is a model parameter whose value is set 
to 1 in the following [5]. Then, the time lag be- 
tween two photons with energy difference AE and 
emitted at the same time is given by 
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and with £! tot = + £1m = 1, = 0.7 and 
H = Tlkms^MpcT 1 . 

This formalism is different from the one used in 
[2, 3], with the addition of an extra factor (1 + z) 
in the integral of Eq. 1 . 

Gamma-Ray Bursts (GRBs) are well suited for this 
kind of study. They are very bright sources at cos- 
mological distances, and their light curves have 
spikes that can be used to measure time lags. In or- 
der to probe the energy dependence of the velocity 
of light induced by Quantum Gravity, we analyse 
the time lags of 15 GRBs observed by the HETE-2 
satellite [6] as a function of the redshift and look 
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Figure 1 : Light curve of GRB 4 1 6 for the en- 
ergy band 8-30 keV before (top) and after (bottom) 
denoising. 

for a dependance of At as a function of z: 

<At> = &Ki{z)+b(l + z), (2) 

where 




and a and b parameters stand for extrinsic (Quan- 
tum Gravity) and intrinsic effects, respectively. 

As in [2, 3], we use wavelet analysis [7] to study 
the light curves. Unlike Fourier transform, this tool 
is well adapted to the study of non-stationary sig- 
nals, i.e. signals for which the frequency changes 
in time. 

In this paper, we will only discuss maxima of the 
light curves. A study of minima is presented in [8] 
along with a more detailed description of the ana- 
lysis. In the next section, we describe the method 
we use to study the light curves and to measure the 
time lags. Then, we give our results in terms of 
a limit on the quantum gravity energy scale Eqq. 
The intrinsic effects are related to the source itself. 
They are discussed in the last section. 

Method 

The analysis of the 15 GRBs with measured red- 
shift follows the steps described below: 

• determination of the time interval to be stud- 
ied between start and end of burst. It is de- 
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Figure 2: Distribution of the derivatives of the light 
curve for maxima located at to. The contribution of 
fake extrema is negligible for | A//Ai| < 0.2. 

fined by a cut above the background mea- 
sured outside of the burst region, 

• choice of the two energy bands for the time 
lag calculations, later called energy scenario, 
by assigning the individual photons to each 
energy band. 14 different scenarios were 
considered, corresponding to 14 different 
values of the energy gap between low and 
high energy bands, 

• de-noising of the light curves by a Discrete 
Wavelet Transform (DWT, see Fig. 1) and 
pre-selection of data in the studied time in- 
terval for each GRB and each energy band, 

• search for the rapid variations (spikes) in 
the light curves for all energy bands using a 
Continuous Wavelet Transform (CWT). The 
result of this step is a list of maxima candi- 
dates, along with a coefficient characterizing 
their regularity (Lipschitz coefficient a and 
its error 6a), 

• association in pairs of the maxima, which 
fulfill the conditions derived from studies of 
the Lipschitz coefficient. 

Different cuts and selections are performed at dif- 
ferent steps in the analysis. A cut on the derivative 
of the light curve at the position of each local ex- 
tremum is done to ensure that it is not fake (Fig. 2). 
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Figure 3: Distributions of Aa (top) and S(Aa) 
(bottom). Arrows give the positions of the cuts. 

Then, for each pair (indexed by 1 and 2 for low- 
and high-energy band, respectively), the following 
quantities are computed: 

At = t 2 - h 
Aa = \a.2 — a\ I 



5(Aa) = \]ba\ + 8a\. 

and the following selections are applied: 

Aa < 0.4 
5{Aa) < 0.045. 

These cuts are based on the distributions shown in 
Fig. 3 and are valid for all the energy gap scenarios 
we consider. 

As a result, a set of associated pairs is produced for 
each GRB and each energy scenario. The average 
time lag of each GRB, < At >, is then calculated 
and used later in the study of the Quantum Gravity 
model described in the previous section. The evo- 
lution of the time lags as a function of the redshift 
z allows us to constrain the minimal value of the 
Quantum Gravity scale Eqq. 

Results 

The two parameters a and b of Eq. 2 were fitted 
for the 14 energy scenarios. Both parameters were 
found to be strongly correlated, as shown by Fig. 4. 
However, results of the fits do not show any signifi- 
cant deviation of a from at ±3a, so in the follow - 



Figure 4: 95% contours for a and b for the two- 
parameter fit for the 14 scenarios for maxima. The 
box at the bottom left shows the position of contour 
centers. 

ing, we derive the 95% CL limit on the Quantum 
Gravity scale Eqq- 

For this, we define a likelihood function L by 

L = CXP ( V — J' 

where M is the energy scale and where x 2 (M) = 



E ■ 

all GRBs 
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where the index i corresponds to each GRB. The 
dependence of on M as predicted by the consid- 
ered model of Quantum Gravity is given by: 



i(M) = 



1 A < E >i 



M 



where A < E >, is evaluated for each GRB using 
the relation 

A<J5> = <J5>2-<-B>i. 

Assuming a universality of the intrinsic time-lags, 
the average value b (and its error a^) of Eq. 3 was 
obtained as the weighted mean of the values bk 
(and their errors Ufc) from the two-parameter lin- 
ear fit: 

b = and a l = 7^ 
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Figure 5: Evolution of \ 2 as a function of M. 

where the index k corresponds to each scenario and 
Wk = For maxima, we obtain b = 0.0023 ± 

0.0027. 

Figure 5 presents the evolution of x 2 (M)/ndf 
around its minimum for maxima of the light 
curves. All scenarios fulfill the condition 
X^ in /ndf < 2. No significant preference of any 
value of M is observed for most of the scenarios. 

The 95% CL lower limit on the Quantum Grav- 
ity scale was set by requiring x 2 /ndi to vary 
by 3.84 from its minimum. The values of lim- 
its on Eqg were obtained for the 14 energy sce- 
narios. In good agreement with the slope pa- 
rameter a, all values for maxima are within the 
10 14 -10 16 GeV range. The scenario which pro- 
vides the strongest constraint on the QG scale 
(8-20 keV, 60-350 keV) provides us a lower limit 
of E QG > 3.2 x 10 15 GeV. 

Discussion 

GRB light-curves are not "perfect" signals, since 
they exhibit intrinsic time lags between high and 
low energies. The peaks of the emission are shorter 
and arrive earlier at higher energies [9, 10, and ref- 
erences therein]. These intrinsic lags have a sign 
opposite to the sign expected from Lorentz viola- 
tion. The effect of Lorentz violation must therefore 
be searched with a statistical study analyzing the 
average dependence of the lags with z, and not in 
samples limited to one GRB. 



A strong anticorrelation of spectral lags with lumi- 
nosity has been found in [11]. At low redshifts, we 
detect bright and faint GRBs which present a broad 
distribution of instrinsic lags whereas at high red- 
shifts, we detect only bright GRBs with small in- 
trinsic lags. This effect could mimick the effect 
of Lorentz violation due to the non uniform dis- 
tribution in luminosity of the GRBs in our sam- 
ple. Here, with limited statistics, we made a test by 
performing our analysis on a restricted GRB sam- 
ple, almost homogeneous in luminosity. The ob- 
tained results show lower sensitivity (lower values 
of the limits on Quantum Gravity), because of a 
decreased statistical power of the restricted sample 
and a smaller lever arm in the redshift values. 

Conclusion 

Light curves of 15 GRBs with known redshifts ob- 
served by the satellite HETE-2 have been studied 
using wavelet analysis in order to look for a Quan- 
tum Gravity effect on light propagation. No effect 
is detected above ±3cr and a lower limit on Eqg 
is set to 3.2 x 10 15 GeV. This limit can be consid- 
ered as competitive considering the energy gap of 
- 1 30 keV provided by HETE-2. 
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